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Abstract 

We find that spin polarized disordered Fermi liquids are unstable to the nucle- 
ation of superconducting pairing states at mesoscopic scales even in magnetic 
fields substantially higher than the critical one. We study the probability 
of finding superconducting pairing states at mesoscopic scales in this limit. 
We find that the distribution function depends only on the film conductance. 
The typical length scale at which pairing takes place is universal, and de- 
creases when the magnetic field is increased. The number density of these 
states determines the strength of the random exchange interactions between 
mesoscopic pairing states. 
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The stability of a superconducting state with an order parameter A can be characterized 
in term of the generalized curvature in the following way, 

K ' ' 5A*(r)5A(r')' {) 

The curvature determines the stability to the spatial variation of the modulas of the order 
parameter and the stability to the creation of supercurrents. Here E({A(r)}) is the energy 
of a configuration {A(r)}. The curvature evaluated at the ground state should be positive 



defined, i.e. det 0(r,r') > 0. For a dirty superconductor where sjD/A 3> I 3> pp 1 , the 
curvature at A(r) = A has mesoscopic fluctuations, like other physical quantities' 1 '. Here D 
is the diffusion constant and / is the mean free path, p F is the Fermi momentum. However, 
in the absence of an external magnetic field, the mesoscopic fluctuations are small. Thus 
the curvature is almost positive defined, and the conventional homogenous superconducting 
state is stable. 

When the magnetic field is applied parallel to the disordered thin superconducting film, 
the suppression of superconductivity is mainly due to Zeeman splitting of electron spin en- 
ergy levels' 2,3 '. In the strong spin-orbital scattering limit r so A <C 1, close to the critical 
field H®, the average spin polarization energy E p is nearly equal to the average conden- 
sation energy E c . Mesoscopic fluctuations of spin polarization energy due to mesoscopic 
fluctuations of spin sucseptibility become comparable or larger than the energy difference 
E c — E p . Therefore, the amplitude of the mesoscopic fluctuations of the curvature evaluated 
at A = A(H) could be comparable with the average. Here A(H) is the order parameter at 
given H, A is the order parameter for H = and l/r so is the spin-orbit scattering rate. 
At if° — H/H® ~ g~ 2 , where g = e 2 v Dd/(2ir 2 h) is the dimensionless film conductance, u 
the density of states and d the film thickness, the curvature, averaged over the area of size 



£(H) = ^ DAq/ A 2 (H), has a random sign, and the system is unstable with respect to the 
creation of normal regions or spontaneous creation of supercurrents' 4 ' . 

One of the consequences of the mechanism discussed above is the instability of the spin 
polarized disordered Fermi liquid well above the critical magnetic field. This is because 



though the average curvature of the normal metal state (0(r, r') evaluated at A = 0)is 
positive defined, i.e. E p > E c , its mesoscopic fluctuations have random signs because of 
the mesoscopic fluctuations of the spin polarization energy. In the regions where the spin 
polarization energy cost to form superconducting pairing state is much lower than the average 
energy cost, the fluctuations of the curvature are of large negative value comparable to its 
positive average such that the normal metal with A = becomes unstable. As a result, above 
the critical field H®, the superconducting pairing correlations are established at mesoscopic 
scales in the different regions in the normal metal and couple with each other via exchange 
interactions of random signs. 

In this paper, we study the probability to find regions where the superconducting pairing 
states are formed at mesoscopic scales at H > if°. At high magnetic fields in the strong 
spin-orbit scattering limit, the statistics of these pairing states can be studied with the help 
of the generalized Landau-Ginsburg equation, 



e (V - i|A) + A(r) + / 50(r, r', H) A(r')dr' = (2) 

where £ = \J ' D/ A . The statistical property of the random potential SO(r, r) is determined 
by its second moment in the Gaussian approximation M 



(SO{T U r[)SO(T 2 y 2 )) \ 



£o 2 5( ri - r 2 )5(r; - r' 2 )5( ri - ri) + 5(n - r^ft - r 2 ) 



(3) 



where (...) denotes the average over the impurities realizations. The curvature evaluated at 
the normal metal state where A(r) = has a simple form: 



0(r,r>) = - 



, . 2e \ 2 H°-H 



5(r-r') + 60(t,t') (4) 



c J H° 

Eq.2 is a nonlinear equation in terms of A(r), with a nonlocal SO(r,r') potential orig- 
inating from the oscillations of the wave functions of cooper pairs. Generally speaking, 
it is qualitatively different from the Schroedinger equation of an electron in the presence 
of random impurity potentials^ 5-8 !. These complications arise naturally in the study of 



the interplay between the mesoscopic effects and the superconductivity and are the generic 
features of strongly correlated mesoscopic systems. In fact, this nonlocal structure of the 
potential in Eq.2 leads to the superconducting glass state. 

However, at H — H® 3> H®/g 2 , the optimal configurations which determine the macro- 
scopic properties of the sample turn out to be the superconducting droplets embedded inside 
the disordered Fermi liquid, with the phases of each droplet coupled via random exchange 
interaction. Such a configuration can be characterized by three parameters: A), the typical 
size of the droplet, L f ; B). the typical distance between the droplets, L d 3> Lf\ C). the 
typical value of the order parameter inside each droplet. In the following, we will discuss 
the statistical of the mesoscopic pairing states in this regime. In this limit, in the leading 
order of (Lf/L d ) 2 the statistical property of the formation of superconducting pairing states 
at mesoscopic scales is similar to that of the impurity band tails' 5 ~ 8 l 

Eq.2 has a nonzero solution only in the regions where the fluctuation of the curva- 
ture is of order of average, i.e. H — H®/H® and the corresponding curvature at A = 
is negative. The nonlinear term in Eq.2 determines the amplitude of A and its effect 
will be discussed only when the typical amplitude of A is concerned. The most prob- 
able configurations are those with superconducting droplets embedded inside the normal 
metal, i.e. A(r) = J2a ^aVa{ r ) > / dvq a {r)r]p{v) ^ <W- Note 77 (r) introduced in this way 
is dimensionless. For such a configuration to have lower energy than the normal state, 
/ drdr'A(r)0(r, r')A(r') < 0. The total energy of such a configuration consists of terms cor- 
responding to the coupling between different droplets. The coupling between the droplets 
decays as distance increases. When the size of the droplets is much smaller than the dis- 
tance between them, the typical magnitude of the coupling between different droplets is 
much smaller than that of the coupling within one droplet, provided L d ^> Lj. We are going 
to neglect such terms in the estimate of the probability of the droplets in the leading order 
of o(L 2 /L d ). Thus, to have / droplets in the normal metal, / independent inequalities have 
to be satisfied 
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J dr Va (r) ^oV 2 + ^T^j Va(r) + J ^rV(r)^(r,r')^(r') 



<0 (5) 



Furthermore, we can write down the probability to have superconducting pairing states 
at H ^> H® in term of the sum of probability to have certain number of droplets 
V({r](x)}) = J2i Pi{{f]a}\ot = To simplify the notation, we introduce Olg = 

QV 2 + H- H%/H%, K M = SO(r, r'). Taking into account Drj(r) = U a Dr] a , we have 

Pi({Va}\a = 1, l)=f P({K M })U a N l J 9(-L a + F a )D Va DK M (6) 

where L a ({r](r)}) = J drr] a (r)0 LG r] a (r), F a ({r](r)}) = J drdr f r](r) a K M (r, r')r) a (r'), and TV is 
a normalization constant. We use the following equality to transform the step function into 
integrals: 

/0 z'+oo 
dg a / dh a exp [ih a (L a -F a - g a )\ . (7) 
-OO J —CO 

Eq. 6 is reduced to Pi({rj a }\a = = Yi a p a . In the Gaussian approximation, when 

the statistics of 50(r, r') is completely determined by the second moment of the correlation 
function, p a can be simplified in closed form as 

J drr] a (r)0 LG r] a (r) 



p a = N J erfc 



^fdr^ j dr 2 dr' 2 C{r 1 , r[; r 2 , r / 2 )r/ a (r 1 )r7 Q (r / 1 )77 a (r 2 )^(r / 2 ) 



Dr) a (8) 



where C(r 1: r[; r 2 , r' 2 ) = (5050) as given in Eq.3. erfc(a/6) = J a exp(— x 2 /2b 2 ) / V 2nb 2 dx. 
One can evaluate the functional integral Dr) a (r) in the saddle point approximation as long 
as H — H® 3> H®j 'g 2 . The saddle point equation of Eq.8 can be obtained by minimizing the 
argument of the error function. The solution of the saddle point equation i] s (r) determines 
the shape of the optimal droplets. To carry out the functional integral of r) a (r), one can 
expand rj(r) around the saddle point, rj(r) = rj s (r) + 5r](r), Srj(r) = J2 n ^(r), where ry n (r) 
are the eigenstates of the operator T(r, r') generated via second functional derivative of 
the argument in the error function with respect to ^(r) at 77(1") = ^ s (r). Our final result 
depends slightly on the detail structure of T(r, r') and we do not give an explicit form here. 
Performing the gaussian integral of 5rj(r) around the saddle point, taking into account the 
normalization condition, we obtain, 



„ . L s . detTfr, r') 
p a = erfc(-=g=) — — ) / [doo] (9) 



fWJ det (C(r,r')) 

where L s / v / 2a s is the argument of error function in Eq.8 evaluated at 77 (r) = rj s (r). ' 
indicates the exclusion of the zero eigenvalue. The last integral in Eq.9 corresponds to the 
contribution from the zero eigenvalue state, originating from the translation invariance of 



the saddle point equation, with 2 — fold degeneracy r] 0i (r) = L d i7] s (r — r )/yj r] 2 dr,i = 
x, y' 6 ' 7 l Here L is the characteristic length of the droplets determined via the normalization 
condition 1 = / drL^V^) 2 / 2 IVs dr - Th us, 

/ [da ] = -L I dx a dy a . (10) 

J -Lq Jv a 

The spatial integral is performed only in the region v a where no other droplets are present. 
Using the following rescaling 

r = yL f , V = V y Lj\ r]s (r)=r ]s (-^), C(r, r'; n, r[) = ^C(y, y'; Vl , y[), L f = £ ' 



we can express L s , a s in term of dimensionless r) s (y) 

T - p H ~ H c r 2 rr - A 2 ^l 
H» 1 g 2 

B = J dyr] s (y)(V 2 + l)r] s (y), A 2 = J dy 1 dy[dy 2 dy' 2 C(y 1 ,y[]y 2 ,y f 2 )r] s (y 1 )r] s (y[)r] s (y' 2 )r] s (y 2 ) (12) 

where B, A 2 are the dimensionless quantities of order of unity depending on the details of 
rj s (y). r] s satisfies the dimensionless saddle point equation 

(Vj + l)Vs(y) + / dy 1 dy[dy'C(y,y';y 1 ,y[)r ]s (y 1 )r ]s (y[)r ]s (y') = (13) 

and at y = 00, rj s (y) = 0. We also estimate that det' T(r, r')/det (0(r, r')) ~ I. 
Collecting all the results, we have 



Pi({Va}\a = oc — |^2 erfc 



V 1 f 1 „ \Bg (H C -H\ 1/2 Y' 



H, 



c 



(14) 



where V l /l\ is from the spatial integral in Eq.10, excluding the overlap between different 
droplets. We take into account L ~ Lf. It is easy to confirm that the average number 
density of the droplets is 
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Bg ( H-H* \ 1/2 ' 

A \ H c J 



(15) 



Let us turn to the problem of the typical amplitude of A a of a droplet. The amplitude 
of A Q is determined by the nonlinear term of Eq.2 and the probability to have a supercon- 
ducting droplet with A Q = A is 

2AiV r / A \ 2 

V(A) = J P(K M )5 N a {—) +L a -F a dK M d Va (16) 

where N a is given as N a = J rj a (r)dr. L a , F a are given after Eq.6. The prefactor in front of 
the integral is from the Jacobian under the transformation S(x) — > 2AxS(Ax 2 ). Transforming 
5-function into a integral and carrying out the gaussian integral of Km , using the procedure 
similar to that used in obtaining Eq.8, we can expressed Eq.16 in terms of rj a (r). The 
functional integral can be performed in the saddle point approximation. This saddle point 
equation turns out to be the same as that of Eq.8 except there is an additional nonlinear 
term proportional to N a (A 2 /Al). However, as shown below, the typical A in the optimal 



droplet is much smaller than Aoy H — H®/H® in the limit H — H®/H® ^> 1/g 2 . Therefore, 
this nonlinear term is much smaller than the linear term already present in Eq.13 and can 
be treated as a perturbation as far as the spatial dependence is concerned. As a result, we 
can use the saddle point solution obtained in Eq.13 to evaluate Eq.16. By expanding the 
resultant equation in term of the nonlinear term (A/Ao) 2 and keeping only the A dependent 
term, we obtain the conditional distribution function of A of a droplet 



^ c (A) = -^exp(-^ 2 TI ). (17) 



2CVA / 2 A 2 ' 
^ r exp 

This shows that A ~ A /g is independent of the magnetic field. C is a constant of order of 



unity. It is much smaller than A y H — H®/H®, justifying the approximation we made to 
derive Eq.17. 

So far we neglect the coupling between different droplets and treat the droplets as a 
dilute gas. Though the coupling between droplets does not affect the probability of finding 
one droplet, it determines the coupling between the different droplets and the global phase 



rigidity. The typical distance L d is order of L f eiic~ l/2 (Bg/A(H - # c °/# c °) 1/2 ) following 
Eq.15. The typical coupling between two droplets is determined by SO and is given as 

z/ A Q A /3 J drdr'0(r, r')r] s (r - r a )r] s (r' - r p ) 

To obtain this result, we take into account that the size of the droplet is Lf, typical A a 
is given by Eq.17, and |r a — rg| ~ L d . The average O, as shown in Eq.3 is proportional 
to 5(r — r'). Thus, the average coupling is proportional to the overlap integral / drrj s (r — 
r a)Vs( r ~ Y p) ex P(~Ld/Lf), which is small in the limit L d ^> Lf. This indicates that 
the distribution function of the coupling between different pairing states is symmetric with 
respect to zero. 

The existence of random Josephson coupling in the presence of a parallel magnetic field 
is a consequence of the Pauli spin polarization. This phenomena exists even without spin 
orbit scattering. Consider for example a granular superconductor, with superconducting 
grains coupled with each other via Josephson coupling. The sign of the Josephson coupling 
is determined by the total phase of the time reversal pairs. In the pure limit, though the sign 
of the wave function of each electron oscillates with a period of Fermi wave length, the total 
phase of (p, — p) pair is zero because of the exact cancellations of the phases of each electron 
inside the pair. Therefore there is no sign oscillation for Josephson couplings. In the dirty 
case p is not a good quantum number. However the sign of the coupling doesn't oscillate 
as a function of spatial coordinate because of the time reversal symmetry. As a result, even 
when the distance between the grains is much larger than the mean free path, the sign of 
the coupling is positive definite' 9 !. This is in contrast to the ferromagnetic ordering nuclear 
spins due to RKKY exchange interaction. In the pure limit, RKKY coupling exhibits Friedel 
oscillations with the period of Fermi wave length. In the presence of impurity scatterings, 
the phase of Friedel oscillations of electron wave functions becomes random. As a result, 
the nuclear spin system exhibits spin-glass type ordering instead of ferromagnetic ordering 
when the impurities are present. 

In the presence of a parallel magnetic field, the electrons inside the normal metal become 
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oc — -erfc 
Q 




(18) 



polarized. In this case, the electron with spin up has a different energy as the electron with 
spin down at the Fermi surface because of the Pauli spin polarization. As a result, the phase 
of the electron spin up does not cancel with that of spin down in the presence of Zeeman 
splitting, and the total phase is equal to J p up • dr + J pdown • dr — J drfi B H/vF- Here 
Hb is Bohr magnetor and v F is the Fermi velocity. This leads to the sign oscillations of the 
Josephson coupling. We assume the electrons in the metal are fully polarized but neglect the 
spin polarization effect inside the grain because the Zeeman splitting energy scale is much 
smaller than the energy gap inside the grain. More specifically, Josephson coupling between 
these two superconducting grains is proportional to the cooper pair correlation function^ 10,11 ! 
TJ2 n { a y^G^(r 1 , r 2 )a^ 5 G s " 6n (r 1 , r 2 )}, where a y is the y component Pauli matrix, Gfj is the 
Green function in Mastubara representation with spin index /3, 7. In the clean limit when 
I ^> ^bH/vf at T = 0, the cooper pair correlation function is 

[^n d2 Q cos [Q • (n - r 2 )] . e ( 2fi B H tt\ 1 

/ dedO- — — - 1 — tanh — — ) oc cos r x — r 2 - - 19 

J (2n) 2 2e + 2fi B H + v F • Q K 2kT J V 1 v F 4y|n-r 2 | 2 v ; 

for 2D case at large distances. 9 is the angle between Fermi velocity and Q. The sign of 
the coupling oscillates with a period vf/^bH, which is much longer than the Fermi wave 
length, with which the sign of RKKY interaction oscillates. 

In disordered metals, for electron pairs to travel between grains separated with a dis- 
tance L, they must typically move along paths of length of order of L 2 /l. When L 2 /l <C 



vf/^bH, the sign is unpredictable. Indeed in the dirty limit when L ^> \J D/ fi B H I, 
the Josephson coupling averaged over the impurity configuration is exponentially small 
exp(— \/2\r 1 — r 2 \/ \J D / n B H) while the typical amplitude of the coupling decays in the 
same way as in the pure limit, i.e. (\Ej\) oc |ri — r 2 |~ 2 . Therefore when the magnetic 
field increases, only the position of the maximun of the distribution function moves towards 
zero while the width of the distribution function barely changes. This results in supercon- 
ducting glass states. Note that in principal the charging effect inside the grain will also 
lead to superconducting glass phase as suggested in a recent experiment ' 12 '. However in the 
metallic limit when the tunneling conductance between the grain and the normal metal is 



much larger than e 2 /h such an effect is negligible. 

It is worth pointing out the mesoscopic pairing states we discussed in this paper are not 
originating from the inhomogeneity of the impurity concentrations. In fact, the range of 
the impurity potential is of atomic scale due to the perfect screening in the metal and the 
fluctuation of the impurity concentration within the area of coherence length is inversely 
proportional to \fn im d£ 2 , and becomes negligibly small in the disordered metal. Here rii m 
is the impurity concentration. Most importantly, the long range coupling mentioned above, 
which leads to the superconducting glass state, is not related to fluctuations of any local 
quantities. We are very grateful to B. Altshuler, L. I. Glazman, D. Huse, I. E. Smolyarenko, 
B. Spivak for useful discussions. F .Zhou is supported by Princeton University. C. Biagini 
is supported under PRA97-QTMD in Italy. We also like to thank NEC Research Institute 
for its hospitality. 
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